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The effects of meson fluctuations are studied in a nonlocal generalization of the Nambu-Jona- 
Lasinio model, by including terms of next-to-leading order (NLO) in 1/Nc. In the model with only 
scalar and pseudoscalar interactions NLO contributions to the quark condensate are found to be very 
small. This is a result of cancellation between virtual mesons and Fock terms, which occurs for the 
parameter sets of most interest. In the quark self-energy, similar cancellations arise in the tadpole 
diagrams, although not in other NLO pieces which contribute at the ~ 25% level. The effects on 
pion properties are also found to be small. NLO contributions from real mv intermediate states 
increase the sigma meson mass by ~ 30%. In an extended model with vector and axial interactions, 
, there are indications that NLO effects could be larger. 
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Until recently most calculations performed with the Nambu-Jona-Lasinio (NJL) model [y^] have been restricted 
to tree-level in the mesons or, in other words, to leading order (LO) in the 1/Nc expansion. Although the possible 
] importance of quantum fluctuations of the mesons has been generally recognized (see, for example, Refs. |3|-|^), a 
satisfactory treatment of such effects has required the development of schemes that are faithful to the chiral invariance 
of the model. This is crucial, given that the most important motivation for studying the model is its description of 
the dynamical symmetry breakdown. Schemes have been introduced on the basis of an effective action method 
' an appropriate selection of Feynman diagrams and a bosonized approximation ||] . 

, Calculations at next-to-leading order (NLO) in l/N^ are certainly more involved than at LO, both analytically 
' and numerically. However, since 1/Nc is such a modest expansion parameter, it has rightly been seen as important 
to estimate the size of some NLO effects. Even if the calculations were of interest for no other reason, a check on 
this perturbative approach would be very valuable. Encouragingly, there seems to be a general consensus emerging 
pL( from recent studies that NLO corrections are indeed relatively small for standard choices of the model parameters 
IQ-pq], at least for properties of the ground state and of the pion. A contrary view, however, has been expressed 



I. INTRODUCTION 



I 

by Kleinert and Van den Bossche |1^ . 

(~| , The most compelling reason for a NLO analysis, however, is that the LO treatment of an NJL-type model neglects 

^ ■ physical processes that are expected to be qualitatively important. For instance, several of the particles described by 

• ^ , such models (tr, p, ai) are broad resonances, yet the model meson propagators at LO are purely real. At NLO the 

' particle widths are incorporated in a natural way, by including diagrams with two-meson intermediate states in the 

H " Bethe-Salpeter equation (BSE). Although these widths can be estimated from the LO decay amplitudes, the effects 

. . 1 on the position of a resonance require a NLO treatment. Such effects may well prove important in model descriptions 



of, for example, the p meson |}4|,|13[. 

A difficulty with the model at NLO concerns the specification of model parameters associated with regularization 
of the loop integrals. At LO the model contains three parameters: the current quark mass, a coupling constant 
and a cutoff parameter that regularizes quark loops. Conventionally two of these parameters are fitted to the pion 
mass and decay constant, while the third is at least constrained by the value of the quark condensate. However, 
the regularization of the model must be further specified at NLO since meson and quark loop integrals can be 
regulated quite independently without affecting chiral symmetry relations [|6|J^j9|,|l^- p^ , p7t . The parameter space can 
be limited to some extent because the NLO corrections produce instabilities for large values of the meson cutoff p3| , |l7| . 
Nonetheless the freedom to introduce another cutoff parameter produces significant uncertainty. Oertel et al. ||13| ] 
have determined the meson cutoff by studying the pion electromagnetic form factor and find values well away from the 
region of instabilities. This is a promising result, but to put it on a firmer basis one would have to relax assumptions 
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that were made in Ref. ||T^ for numerical reasons^. An attractive alternative is to perform the regularization at the 
level of the model action rather than in the loop integrals. In this approach the regularization is specified from the 
outset and does not have to be imposed at each order of the expansion, a point which has been stressed by Ripka Q . 

In this paper, we make a NLO analysis of a model that was originally proposed in Ref. It has subsequently been 
studied at LO in the meson baryon (soliton) ||2^ and quark matter sectors. The two-body interaction 

vertex is nonlocal and separable, a Gaussian form factor being associated with each quark field. Such a model 
successfully eliminates several of the traditional problems of the NJL model whilst nevertheless retaining much of the 
simplicity that is its chief merit. The separable nature of the interaction is motivated in part by instanton liquid 
studies ||2^,^. The interaction form factors ensure the convergence of all loop integrals and so the NLO corrections 
are unambiguous. In addition, a practical advantage is conferred by the shape of the form factor, which allows 
complicated NLO diagrams to be evaluated efficiently with Gaussian numerical techniques. Such properties make this 
model a particularly convenient one with which to examine NLO effects. NLO quantum fluctuations of the quark 
condensate in a variant of this model were studied in Ref. Here we investigate the NLO treatment of mesons. 

In our approach we keep the full momentum dependence of all quark loop diagrams, and avoid expanding them in 
powers of meson momenta. Such an expansion may be useful in the determination of pion properties, but it is not 
valid for heavy mesons, such as the scalar, isoscalar sigma. 

II. THE NONLOCAL MODEL 

In this section we define the model used and recall some important aspects of its behaviour at LO. This is entirely 
standard Pp9|] and is included largely in order to establish our notation. The action of a quark model with a four-point 
interaction may be written 

s= /A?,.)w-,„d,H.)+i: /nA,.H.(.„.„.3,.,) 

X i^{xi)Tfi:{x3Mx2)T,^ij{x4), (1) 

where Tf denotes Dirac, colour and isospin matrices. Chiral symmetry constrains certain of the possible Dirac and 
isospin structures to appear in particular combinations. The admissible interaction terms are 

HsiY' 7m). ^^4(7''75 » 7^75), 
i^5(r"0r'^+^75«'^75), iJsK. ® - 't^.t'^ ® ct^'t'^). (2) 

We shall be primarily concerned with the simplest version of the model, including the pions, and their chiral partner, 
through the Hi interaction. In the usual NJL model the Hi are simply constants. However, the use of a nonlocal 



interaction is appealing [g^ 19 2l|], not least because it enables one to avoid some of the ambiguities associated with 
regularization. Moreover, a nonlocal interaction of separable form is suggested by studies of the instanton-liquid 
picture of the QCD vacuum [p2|,p3|. In this case, the momentum space form of the Hi is 



Hi{pi,P2,P3,P4) = 5(27r)''G,/(pi)/(p2)/(P3)/(P4)<5(pi +P2 -P3 'Pa), (3) 

with the Gi being constants which, for the purposes of counting, we take to be 0{1/Nc). Although the instanton-liquid 
model (2^:^ predicts a particular form for the form- factor function f{p), other authors ||l8|,|l9|,^ have considered a 
function which is a Gaussian in Euclidean space, 

fiPE) = exp(-p|/A2). (4) 



The Gaussian shape produces a good phenomenology of mesons and nucleons in the LO approximation [ [l9| , pi[ and is 
used for the numerical results in this paper. 



'^The p meson was included in the model but vector and axial- vector degrees of freedom were not included as intermediate 
states in the NLO diagrams. In addition, nai mixing was neglected, even though this is a LO effect. 
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At leading order in the 1 /Nc expansion the quark propagator is dressed by a single quark- loop self-energy diagram. 
This results in a momentum-dependent quark "mass", m{p): 

S-\p) = i> - m{p)- (5) 



m[p) = mc+ {rn{Q) - m^f[p), (6) 

where the constant to(0) is obtained from the Schwinger-Dyson equation (SDE) [ p8[ . If m(0) is sufficiently large, the 
poles in the quark propagator are shifted into the complex p^ plane A feature of models with nonlocal 

interactions is the appearance of additional, unphysical poles in the quark propagator. For the parameter sets of most 
interest these are located well away from the real axis [|9| , but instabilities of the model ground state can appear for 
small values of the cutoff |Q. 

Mesonic bound states can be found from the poles in the qq scattering matrix, T. For a separable interaction, it is 
convenient to factor out interaction form factors of the external quark momenta and to define 

T{pi,P2,P3,P4) = S{pi +P2 -ps - P4)Y[ f {P7l)^^iT^J{q)T j , (7) 

where q — pi — ps — Pi—p2 denotes the total momentum of the qq pair. The notation Tj is needed in the longitudinal 
components of the vector and axial channels where, for example, Tv = — Ty = i^. In all other cases T j and Tj are 
the same (see Eq. (23) of Rcf. ^). The Bethe-Salpeter equation (BSE) at LO is given by 

f{q) = G + GJ{q)f{q), (8) 

where G is simply a matrix of the coupling constants from the action and J{q) is the matrix of polarization loop 
integrals 

J,, (g) = zTr j ^f{p^)f{p^)T,S{p^)T,S{p+), (9) 

where we have introduced p± = p ± ^q. The symbol 'Tr' is used to denote a trace over flavour, colour and Dirac 
indices. 

Close to the pole corresponding to a particular particle, the amplitude in the relevant channel may be written 

,10) 

m — q 

where V{q) denotes the particle vertex function. In the simplest form of the model, with only the Gi coupling 
included, these functions are: 

V^{q) ^ g^qqil^r'' , Va{q)=gaqq, (H) 

where: 



1 _ dJii 
9%q ~ dq^ 



(12) 



The vertex functions appear in calculations of mesonic properties, such as the pion decay constant, /^r. 

In order to obtain symmetry currents with the same divergences as the corresponding local currents in QCD, and 
hence to respect the corresponding Ward identities, one has to include additional nonlocal terms in the currents 
| p^ , p^ 21 2^. These arise as a consequence of the nonlocality of the action. Details of the nonlocal terms required in 
both vector and axial currents can be found in Refs. |l^,|l^. Although the method developed in those papers is not 
unique and there is an ambiguity in defining the transverse parts of the currents, this does not affect the longitudinal 
components. 

As an example, consider the pion decay constant. The diagrams in Fig. ^ represent the contributions from the local 
and nonlocal pieces of the axial current. The longitudinal part of the current, contains the following term 
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2(2^)12 J Y{'^'^Pni^iPl)h5T°'i^{p3)lp{P2)^{P'l) S{pi + P2 + q - P3 - Pi) 

X [f{pi)f{p2)f{p3)f{p4 -q) + f{pi)fiP2)f{p3 - q)f{Pi) 

-f{pi)f{p2 + q)f{p3)f{p4) - f{pi + q)f{p2)f{p3)f{Pi)] ■ (13) 

This term contributes to the pion decay constant since the operator V'(P2)'*/'(P4) has a nonzero vacuum expectation 
value (represented by the bubble loop in the second diagram of Fig. p. The rest of the operator, ipiPi)il5T°"4'iP3)j 
contributes to the pion-to- vacuum matrix element of the current. Such terms are needed in order to satisfy the 
Gell-Mann-Oakes-Renner (GMOR) relation ||2^, as discussed in Ref. and they can make significant numerical 
contributions to various observables ||l9[| . 



III. NLO CONTRIBUTIONS 



A. Symmetry currents 

At leading order in the l/N^ expansion the fields at Xi and 0:3 in the interaction term of the action (|l|) are always 
connected to the same quark loop, and similarly for the fields at X2 and X4. Beyond this order, there are also exchange 
or "Fock" terms where the fields at 0:1 and X4 are connected to the same loop. These can be constructed by first 
performing a Ficrz transformation on the action and then using the transformed action in just the same way that one 
uses the original action at LO. The Fierz transformed action of this model contains the following colour-singlet terms: 

{4Nc)-\Gi - 2G3 + 2G4 - G5 + 12G6) (10 1 + ^75T" ® ^75T"), 
{m,)-\-2G2 +G3 + G4) (TpT'^ ® 7^r'^ + 7^75t" 7^75t'^), 
{m,r\-2Gi + 6G2 + Ga + G4 - 2G5) (7m ® 7"), 
(4iVe)-'(2Gi + 6G2 + G3 + G4 + 2G5) (7A.75 «'7^75), 
i4N,)-\~Gi - 2G3 + 2G4 + G5 - 12G6) (r" ® t'^ + ^75 ^5), 

{8N,)-\Gi - G5 - 4G6) (a^. ® a'"- - a^,r" ® a^'^r"). (14) 

(There are also colour-octet terms which we do not use here.) 

The nonlocal terms in the symmetry currents of the model also involve four quark fields and so will also be subject 
to exchange effects at NLO. These Fock pieces in the currents will introduce further ambiguity through the definition 
of their transverse parts. One way to construct the NLO current terms is to take the Fierz transformed action 
and to apply the same method as was used in Refs. [p^p^ to form the nonlocal currents from the original model 
action. This method leads to nonlocal NLO current terms of the same structures as those presented in Ref. ||l^, with 
the appropriate coupling constants in the original terms replaced by the corresponding combinations of couplings 
from the Fierz transformed action. For example, as well as a term Gi(«75T° ® 1) in the nonlocal axial current 
constructed from the original action (see Eq. (|3|)), there is a Fock term of the same structure where Gi is replaced 
by (4Ar,)-i(Gi - 2G3 + 2G4 - G5 + 12G6). 

An alternative approach to constructing the Fock terms of the model's currents is simply to Fierz transform the 
LO currents of Ref. [|l^. The differences between this procedure and the one outlined above (where the Noether- 
like method for deriving a symmetry current is applied after the Fierz transformation) lie in the combinations of 
form factors which are present in the currents. However, these affect only the transverse parts of the currents; both 
definitions satisfy the appropriate Ward identities. The equivalence of the longitudinal components has been checked 
explicitly p^ ]. In practice, the currents obtained from the Fierz-transformed action are somewhat easier to work with, 
since they retain the same structures as the nonlocal terms in the LO currents, with appropriate substitutions for the 
overall coefficients. 



B. Quark propagator 

We are now in a position to consider the quark and meson propagators at NLO. Like Oertel et al. ||l^,|l^,|l^ and 
Ripka 1 14 1 we use a strict NLO scheme in which all propagators and matrix elements are expanded to NLO in l/Nc- 
This scheme has also been presented by Dmitrasinovic et al. ||], who themselves preferred to adopt an alternative 
scheme where certain NLO terms in the quark self-energy are treated self-consistently in the quark Schwinger-Dyson 
equation. Their scheme resums a subset of terms to all orders in l/Nc- Both schemes satisfy chiral low-energy 
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theorems [^,|3[. The Feynman diagrams required by the strict NLO scheme in the nonlocal NJL model are briefly 
presented below. 

At NLO the quark Schwinger-Dyson equation is coupled to the meson Bethe-Salpeter equation and takes the form 

§ 



f{p-k), (15) 



where Tij(k) denotes the element of the qq scattering matrix which describes scattering from the state with matrix 
structure Tj to the state with structure Ti . This matrix is diagonal apart from the mixing between the pseudoscalar 
and axial channels which is present in the extended model with vector and axial interactions. It is to be understood 
as the LO scattering matrix, which is of order (see Eq. (||)). 

Although the SDE as written above contains all of the required terms at LO and NLO, it also includes some 
unwanted higher order terms. The NLO terms can be obtained by replacing the quark propagator in the final integral 
by its LO form, and by expanding the propagator in the quark "bubble" integral to NLO: 

Sf{p) c s{p) + s{p)j:^'\p)s{p), (16) 

where S^^-' denotes the NLO contribution to the quark self-energy. This leads to an equation for S^^' of the form 

I](i)(p) = ^G,f{p)Trj ^^S{k)Y.^^\k)S{k)f{k) 

-^f{p) E / 7Sl4 (fc)r.5(p - k)T,f{p - k). (17) 

The first term comes from an insertion of the NLO self-energy into the quark bubble diagram. The second, involving 
the scattering matrix, includes both Fock terms and dressing of the quark by virtual mesons. 

The separable nature of the interaction means that this equation can be solved straightforwardly to get 



{p-k), (18) 



^,3 

where 



l-GiJ,,(0)^ J (2^)4 y (2^)4 

xSik)T,S{e)TjSik)f\k)fi£). (19) 

The two terms correspond to the tadpole and meson cloud diagrams shown in Fig. |^. The tadpole diagram is 
responsible for the contribution f^(p)C. It includes the exchange of a zero-momentum sigma meson between the 
quark and a virtual meson. Since C is a momentum-independent constant, this contribution has the same form as 
the LO quark self-energy. The other diagram describes the emission and subsequent reabsorption of a virtual meson. 
It has a nontrivial dependence on the momentum of the quark and generates a wave function renormalization as 
well as a scalar term. Note that although wc shall often refer to the upper double line in such diagrams as a meson 
propagator, it is really a scattering amplitude and so the diagrams contain Fock terms as well as virtual meson 
contributions. 

C. Meson propagators 

At LO, the mesonic bound states are constructed from the ladder Bethe-Salpeter equation (Eq. (||)). The NLO 
extension of the BSE can be expressed in terms of corrections to the basic two-quark loopj^of Eq. (^. If the quantities 



^We refer to a loop integral containing n LO quark propagators as an n-quark loop. 
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Jij are redefined to include these NLO contributions, the scattering matrix can still be written in the form of Eqs. (0) 
and d). 

An obvious set of NLO terms in a polarization loop integral Jij (q) consists of insertions of the NLO quark self- 
energies of Fig. 1^ on either the quark or antiquark line. These contribute 



»Tr / j:^T,Sip^)T,Sip+)J:^'Hp+)Sip+)fiP-)f{p+). (20) 



(2^ 

A second kind of NLO contribution arises from the exchange of a i-channel virtual meson between the quark and 
the antiquark. The corresponding diagram is shown in Fig. || and makes the following contribution to Jij{q): 

■^TUp~k)T,S{k^)TiS{p^)T,S{p+)T,^S{k+) 

l,7n 

^fik^)f{k+)f{p-)fip+). (21) 

Finally, there are NLO contributions that involve intermediate two-meson states, represented by the second diagram 
in Fig. ^ These allow for the instability of a meson against two-body decays and so can introduce an imaginary part 
in the corresponding propagator above the threshold energy for the two-particle final states. They are constructed by 
connecting two LO three-meson vertices with two meson propagators. The resulting imaginary part of a meson mass 
can be expressed in terms of the square of the LO decay amplitude. These contributions also generate a shift in the 
real part of the mass which involves a new type of loop integral. 

In writing explicit expressions for these contributions, it is convenient to introduce functions Lij.fe and Lijk to 
describe the LO three-point vertices. The function Lij^k describes the process k — > ij and is defined as 

f d'^p — — 

Lij,k{qi,q2,q) = iT^r J j^--^rkS{p-)rjS{p - ^{qi - q2))TiS{p+) 



f d^p — — 

«Try j^-^rkS{p-)T,S{p + ^{qi - q2))TjS{p^ 



xf{p-)fip+)fip+Hqi-q2)), (22) 



while Lijk describes the process jk 



— f d'^p — 

Li,jk{q,qi,q2) =iTr J j^-^r,s{p+)TjS{p - ^{qi - q2))rkS{p^) 

xf{p-)f{p+)f{p-Hqi-q2)) 



d^p 

xf{p-)f{p+)f{p+Uqi-q2))- (23) 



+«Tr / j^r,S{p+)TkS{p+^,{q,-q2))r,S{p.) 



The contribution to Jij (q) from the two-meson diagrams can then be written in the form 

- ^ X! / ^^L,^kmiq,P+,-P-)fkiip+)f^ni-P-)Lin,jip+,-p-,q). (24) 



Once the full NLO matrix J*^^^ has been constructed, we can use it to find the corrections to meson masses and 
meson-quark couplings. As an example, consider the pion in the version of the model without vector and axial 
couplings. The pion pole is located at 

^-G^[JilHq') + 4lHq')]^0, (25) 
where superscripts (0) and (1) denote LO and NLO terms. The coupling to quarks (Eq. (|T^)) is given by 
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From these we get the squared pion mass to NLO, 



(26) 



— 



dq^ 



o(0)2j(l)/ (0)2-, 



(27) 



and the pion-quark couphng, 



(0)3 

(1) ^ _ g^qq 



j2 7(0) . t(1) 



(28) 



(0) 2 



D. Meson coupling to a current 



In this section, we consider a NLO determination of the couphng between a meson and an external current, taking 
the pion decay constant as an example. The situation is more complicated in a nonlocal than a local model 
because there are two kinds of contribution to be considered, arising from the local and nonlocal parts of the symmetry 
current (see Sec. |l|). 

The LO diagrams are shown in Fig. |^. At NLO several of the extra contributions can be identified straightforwardly, 
either by inserting the NLO self-energy on any one of the quark lines in these diagrams or by using the NLO 
contribution to the pion vertex function. Other contributions are analogous to those in the BSE at NLO. These 
involve either a t-channel meson exchange in the two-quark loop or intermediate two-meson states. They are shown 
in Fig. |. 

The remaining NLO contributions arise solely from the nonlocal piece of the current. The first of these are Fock 
terms arising from the exchange of two of the quark lines at a nonlocal current vertex, as shown in Fig. ^. Since 
they contain only a single colour trace, as indicated by separating the quark lines associated with each of the ipTip 
factors, this diagram is suppressed by one power of Nc compared with the corresponding LO terms. Although, as 
noted above, there is an ambiguity in defining the transverse components of such Fock terms, this does not affect the 
calculation of the pion decay constant. 

Finally there are two NLO contributions which can be thought of as arising from nonlocal current couplings to 
virtual mesons. They are shown in Fig. ^ The first (a) is somewhat similar to diagram ^(d), except that the coupling 
of the nonlocal current to the virtual mesons involves both the quark and the antiquark instead of a separate quark 
bubble. This is analogous to the two-body diagrams that contribute at LO to several of the electromagnetic amplitudes 
described in Ref. p9| ] . The other diagram (b) can be thought of as coupling the current at a meson-quark vertex. 

Although the NLO contributions to the pion decay constant involve rather a large number of diagrams, in practice 
the situation can be significantly simplified as a result of cancellations amongst them |26j . A complete description of 
these cancellations would be somewhat lengthy and so we offer a brief summary in Appendix We have also checked 
the consistency of our treatment by verifying that the GMOR relation |^ holds to NLO (see Appendix ^ . 



IV. NUMERICAL RESULTS 



A. Model parameters 



The multiple integrals involved in the NLO diagrams can be rather time-consuming to evaluate numerically to good 
accuracy. We have therefore chosen to fit the model parameters at LO and then to calculate the NLO changes to 
observables. In fact, for the parameter sets of most interest, we find that the corrections are small for the observables 
used to fix the parameters. Thus, our results should not be very different from those of a full NLO fit. 

We fit the parameters at LO as in Ref. . In the simpler version of the model with no vector and axial couplings, 
we fit two parameters to the pion mass and decay constant. This leaves one free parameter which we characterize 
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in terms of mo(0), the quark self-energy at zero momentum in the chiral Hmit. Details of the sets used are given^ 
in Table |. Values of some quantities calculated at LO with these parameter sets are also given in the table. They 



are qualitatively similar to those obtained in the extended version of the model |19|. For the parameter sets with 
TOo(O) = 300 MeV and larger, the quark propagator has no poles on the real axis and so single quarks do not appear 
as physical states [ p^|JT9| . 

Although most of our work at NLO has used this simpler model, we also show a few results for the extended model 
including vector and axial couplings. For these we use a parameter set with r7io(0) — 300 MeV, which is intermediate 
between the two cases studied at LO in Ref. 

B. Quark condensate 

The NLO term in the quark condensate can be evaluated as the sum of two contributions. These correspond to 
the two self-energy diagrams of Fig. || where the external quark line is formed into a closed loop with a local scalar 
insertion 0. By comparing the values for the condensate at NLO in Table || with the LO values in Table |, we 
see that the NLO contributions are very small. This may seem surprising since virtual pions would be expected to 
make a large (positive) contribution to the condensate. However it should be recalled that the meson "propagators" 
in the NLO diagrams also contain the Fock terms of the two-body interaction. These make a negative contribution 
to the condensate which tends to cancel the contributions of virtual mesons. This agrees with the findings of Ref. 

in a similar model. In the local version of the NJL model, the condensate is also little altered at NLO if proper 
time regularization of the quark loops is used but in the 0(4) scheme there can be appreciable changes at NLO. 
Including only the pole pieces of the meson propagators, significant shifts in the condensate are also found with 0(3) 
regularization |p7| . 

C. Quark propagator 

The quark self-energy at NLO is given by Eqs. (^ and (|l|). Its evaluation requires integrating over the scattering 
matrix T. Evaluating the qq loop integrals J for each point in the mesonic loop integral is very time consuming. We 
have therefore found it convenient to parameterize the J's in terms of a set of smoothly- varying functions. Working in 
Euclidean space, and with the momentum routings of Eqs. (^8|) and (|l9|), the J's always have a spacelike momentum 
argument and so are smoothly-varying functions. We have found that they are well represented by expansions in 
Chebyshev polynomials of the variable x = exp((7^/A^). 

We write the inverse quark propagator as 

Sp\p) = {l + a{p))^-b{p). (29) 

The functions a{p) and b{p) are shown in Figs. and ^, for both spacelike and timelike momenta. They are plotted 
only up to the energy given by the real part of the complex pole in the LO quark propagator. Above that energy 
there are pseudo-threshold effects in the model associated with the continuation from Euclidean to Minkowski space 

At LO a{p) ^ and there is no wave function renormalization. The NLO contributions to the coefficient ^ range 
up to about 0.25, which is consistent with the expected magnitude of 1/Nc corrections. An interesting aspect of the 
results for a{p) is the appearance of a sudden dip just before the pseudo-threshold energy. It would be interesting to 
examine the behaviour of the function above this energy, although that would require a detailed analysis of the quark 
pole structure at NLO, which is beyond the scope of the present work. Also plotted are the individual contributions 
to a{p) which arise from dressing the quark line with virtual pions and with virtual sigma mesons. The pion cloud is 
obviously the dominant contributor. Since its propagator has a pole at small timelike momentum, one expects the T 
matrix in the pseudoscalar channel to be large at modest values of spacelike momenta (the region which dominates 
the NLO integrals). Note also that the pion contributions contain an extra factor of three due to isospin. 

The coefficient of the unit matrix in the quark self-energy at NLO receives contributions from both the tadpole and 
the meson-cloud diagrams of Fig. |^. The results in Fig. ^ show that the tadpole contribution is in fact rather small. 



^Note that the parameters given differ slightly from those quoted in Ref. where a very similar fit was made at the same 
values of mo(0). This is simply because the calculations of Ref. |^ were performed within the chiral expansion of the model. 
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C <C to(0). The meson-cloud contributions are rather more significant, increasing 6(0) by about 25% (again consistent 
with expectations for a l/iVj, effect). The net NLO shift in the quark "mass" at zero momentum &(0)/(l + a(0)) is 
therefore fairly modest, ^ 15%. Quark dressing by meson clouds in the local NJL model has been investigated in Ref. 

It was found that the pion cloud tends to increase b{p) but that this is partially cancelled by the sigma cloud. 
The nonlocal model studied here supports that conclusion and is able to place it on a firmer footing since there are 
no ambiguities associated with the meson loop regularization^ 

Although the tadpole diagrams in the NLO quark self-energy are not numerically significant, C controls the NLO 
changes to the pion mass and decay constant. This can be seen in the chiral expansions of these quantities described 
in Appendix One might therefore wonder whether there is some physical reason why its chiral limit Co turns out 
to be small. For the parameter set with mo (0) = 200 MeV Table || shows that the pion tadpole makes a contribution 
of about 45% of the LO chiral quark mass, but that this is cancelled to a large extent by the sigma tadpole. Such 
cancellations do not persist at larger mo(0), where the pion tadpole changes sign. However for these parameter sets 
both tadpole contributions are small. In these cases it is a cancellation between "true" virtual meson effects and Fock 
terms (similar to what happens in the quark condensate) which is responsible for the small net effect. 

In the alternative scheme of Ref. js), the 1/Nc contribution of the tadpole diagram in Fig. || is resummed to all 
orders when the SDE is solved self-consistently. In contrast, the momentum-dependent contribution of the meson- 
cloud diagram is treated perturbatively. This partial resummation is motivated by an expectation that contributions 
from the tadpole diagrams are larger than the other NLO effects |8|. However we find that such an expectation is 
not fulfilled in the nonlocal model; in fact we find that the contribution of the tadpole diagram is small compared to 
that of the meson-cloud diagram. We thus see no compelling reason to resum this contribution, and so we can avoid 
possible problems with the self-consistent scheme which have been noted in the context of the local NJL model . 

In the extended version of the model there are additional contributions to the NLO quark self-energy from tadpole 
and meson-cloud diagrams with other mesons. Although a calculation of the properties of these mesons at NLO is 
beyond the scope of the present work, it is straightforward to evaluate their contributions to the quark self-energy. 
As in the simpler model, we find that there is only a modest NLO change to the condensate and that the constant 
C is small. Pseudoscalar-axial mixing makes a significant difference to the contribution from the pion cloud to the 
function b{p), but this is largely cancelled by the vector mesons to leave a function very similar to that found in the 
simpler version of the model. 

The function a{p) in the extended model is shown in Fig. ^. Mixing in the pion channel and the clouds of the 
spin-1 states are again important. In a{p) these effects reinforce each other and so the function is significantly larger 
than in the simpler model. An interesting difference between Figs. ^ and || is the absence in the latter case of a dip 
just before the LO pseudo-threshold energy. This is eliminated in the overall result due to the contributions from the 
longitudinal components of the vector and axial propagators. In particular, the steep drop in the pion contribution 
which occurs in the simpler model is removed by mixing with the axial channel. We note that the local NJL model 
has been used in Ref. to describe the p meson at NLO, assuming that the effects of ttoi mixing and intermediate 
spin-1 states could be neglected. The results presented here raise some doubts about those assumptions. Of course 
definitive statements about the quark propagator in the present model cannot be made without a full NLO analysis 
of meson properties in order to refit the model parameters. Since 1/iVc effects could alter meson masses significantly, 
the parameters used here may not constitute a reasonable choice for the extended model at NLO. 

D. Pion properties 

The denominator of the qq scattering amplitude in the pion channel, Eq. (|2^), is shown in Fig. For low energies, 
less than about 250 MeV, the differences between the LO and NLO curves are very small. This implies that NLO 
contributions to both the pion mass and the pion-quark coupling are extremely small. 

Although NLO contributions to the pion mass must cancel in the chiral limit, as described in Appendix such 
cancellations do not influence higher order terms in the chiral expansion. The smallness of the NLO terms in the pion 
propagator is another consequence of the cancellation between virtual mesons and Fock terms, as we saw above in 
the case of the quark propagator. In particular, all of the contributions to the pion mass and decay constant that are 
proportional to the NLO piece of the quark condensate or to C are small. (More details of these quantities are given 
in the appendices.) This supports the usual LO treatment of the pion in this type of model and it also justifies our 



''Note that the tadpole contributions were not included in Ref. Q . Moreover the model meson propagators were approximated 
by the canonical forms for point-like bosons. 
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use of model parameters fitted to the LO pion properties in these NLO computations. Actual determinations of the 
NLO shifts in m^r and /,r have not been made in this work since they are sufficiently small (not more than a couple 
of MeV) that our numerical integration procedures would have to be refined in order to yield accurate values. 

At higher energies, beyond those shown in Fig. ^ NLO contributions in the pion channel do start to become 
significant. When the ttct channel opens, the denominator develops an imaginary part. Shortly before the pseudo- 
threshold energy (twice the real part of the complex energy of the pole in the quark propagator jl^) is reached, the 
NLO terms become sufficiently important to change the qualitative behaviour of the pion Bethe-Salpeter amplitude. 
The slope of the denominator changes sign and there is even a second zero of Eq. (^5|) . This is a potentially worrying 
feature since it corresponds to an unphysical pole in the amplitude, whose residue has the wrong sign. A similar pole 
has been found in the local NJL model |l^. In both models the undesirable high-energy behaviour is caused by the 
insertions of NLO quark self-energies on the quark lines of the basic polarization loop integral. However it is perhaps 
worth noting that this behaviour is reminiscent of that exhibited by a(p) at high energies in Fig. ^ As mentioned in 



Sec. [VC, this behaviour is not present in the model with vector and axial mesons and it is thus possible that the 



unphysical pion pole may be removed by extending the model to include interactions in these channels. 



E. Sigma properties 

The interpretation of the scalar, isoscalar sigma meson in dynamical quark models is complicated by the fact that 
it does not correspond to an experimentally well-determined resonance. The problem, not just for models of such 
a resonance but also for phenomenological determinations of its properties, is its very strong coupling to the two- 
pio n channel. Indeed, while some phenomenological analyses favour a scalar, isoscalar resonance at around 1 GeV 
|p9| , more recently others have found resonances at 600 MeV or lower (for reviews and more complete lists of 
references see Refs. ^l|-^3|). The strong coupling to pions means that any resonance must form a very broad structure. 
Consequently the corresponding pole in the scattering amplitude must lie far from the real axis, which makes its exact 
position hard to determine in a model- independent way |]32| . 

This strong coupling to two pions also means that a bare qq state in the scalar, isoscalar channel should not be 
directly compared with any phenomenological resonance. The effects of the two-pion channel need to be included. In 
the present model, the leading effects of this type arise from the second diagram of Fig. ||. 

At LO in the nonlocal NJL model leads to a qq sigma meson with a mass of less than 500 MeV. This is 

similar to the local model, where the sigma lies on the qq threshold ||l|,|^. Some of the properties of the state are listed 



in Table III. It is strongly coupled to two pions, although its low mass means that its decay width is less than about 
160 MeV. Although wide, this is still narrower than phenomenological fits which tend to give widths of at least 300 
MeV, even for low sigma masses pO|-^. 

The real part of the denominator of the scattering amplitude at NLO is shown in Fig. ^ Above the two-pion 
threshold, there is an imaginary part generated by the contribution to J„„ from the second diagram of Fig. 0. This 
can be estimated from the LO decay amplitude using the Cutkosky cutting rules and it is shown in Fig. |r2| along 
with the results from a direct numerical evaluation. Some care must be taken in the evaluation of this diagram owing 



to singularities of the scattering matrices in the integrand above the decay threshold. We have regulated Eq. (24) by 

replacing T^rn with Gi{l — Giji^J —ie)^^ and linearly extrapolated our results to e = 0, based on several computations 
for e ~ 10-3. 



The values listed in Tabic 111 for the mass of the sigma at NLO are for the mass defined as the zero of the real part 
of the denominator. Although this is a commonly used choice, it might be better to use the real part of the complex 
energy of the pole in the scattering matrix, as discussed by Pennington | |3^ . The choice used here is purely one of 
convenience. In any case the width of the state we find is not so large that we would expect substantial differences 
between the two definitions. 

We find that the sigma mass is increased from its LO value by about 30%. Below the two-pion threshold the largest 
NLO contribution to Jq-o- arises from virtual two-pion intermediate states. This is partially cancelled by the NLO 
quark self-energy insertions to leave a modest net increase in Jo-o- for most parameter sets. On their own these effects 
of virtual mesons would tend to reduce the sigma mass. However, above the threshold, real two-pion states make a 
substantial contribution with the opposite sign, leading to the increase to the sigma mass. This is consistent with the 
qualitative expectation that such states should be important for modelling the sigma. 

NLO corrections to the sigma mass have been calculated by Pallante |^ using a derivative expansion of the bosonized 
NJL model. In that framework the corrections were found to be large and negative, prompting that author to speculate 
that the 1/A^c expansion might break down for the mass of this state. In the local NJL model, the large negative 
corrections to 1 — GiJ^a persist in calculations to all orders in momentum if one uses Pauli-Villars regularization of 
the quark loops and a sharp cutoff for the meson loop momenta ||lq,p4 . The shifts in this case are sufficiently large to 
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move the model sigma pole to negative 9^, producing an unstable vacuum. In both of these calculations, however, the 
results in the scalar sector were very sensitive to the additional cut-off parameter needed to regularize meson loops in 
the local model. In contrast, the magnitude of the NLO mass shift which we obtain in the nonlocal model is consistent 
with expectations for a l/N^ correction. Thus it would appear that the details of the treatment of high-momentum 
states are important in NJL-type models. 

The imaginary part of the polarization loop is shown in Fig. |2| and compared with the estimate obtained by 
applying the Cutkosky cutting rules to graphs with two-pion intermediate states. The difference between these curves 
is caused by the analytic structure of the quark propagator, which contains additional, unphysical poles (Sec. ||). It 
is reassuring to note that these extra poles are located sufficiently far from the low-energy region of interest that they 
do not have a large effect on the imaginary part. For other parameter sets, with large mo(0), the extra poles lie closer 
to the physical region |jl^ and may have some influence on properties of the sigma meson at NLO. Since the model 
predictions may be unreliable in these cases, we have preferred to quote only upper bounds on the sigma mass for 
some parameter sets in Table II], limited by the position of the first pole in the LO quark propagator. 

For energies up to ^600 MeV, the three-quark integral in the (Ttttt vertex function is only weakly dependent on the 
energy. Therefore the variation of Im(Jo-o-) at such energies is primarily a consequence the available two-pion phase 
space. At higher energies the coupling to two pions becomes much weaker^ This is similar to the behaviour of the 
coupling of the sigma to two pions observed in the quark model studied in Ref. p5| (where the scalar mass was taken 
as a free parameter and chosen with the intention of interpreting the model scalar resonance as a heavy, narrow state). 

We note that a two-flavour NJL model of pions and sigma mesons is insufficient for a fully realistic description of the 
spectrum of scalar mesons. One must also take account of strangeness, radial excitations, mixing with qq molecules 
or glueballs, and decays to four pion states. A recent attempt to include such effects can be found in Ref. ||3q|. 



V. CONCLUSIONS 



We have investigated the effects of meson fluctuations on quark and meson properties in a nonlocal NJL model 
[ p^ , p^ by treating it to NLO in l/N^. The two-body interaction between quarks has a separable form, similar to that 
suggested by instanton liquid studies |^,^. At LO the model has been shown to provide successful descriptions 
of mesons and baryons ||l^,^. The interaction form factors regulate both quark and meson loops in a natural way. 
As a result, no new parameters emerge at NLO and the contributions are unambiguous. This makes the model a 
particularly suitable one in which to study NLO effects. 

The contributions to the Schwinger-Dyson and Bethe-Salpeter equations at NLO contain the same structures as 
those in the local model |^,|| , apart from the presence of interaction form factors. The NLO pieces include the dressing 
of quarks and mesons with virtual mesons and the coupling of mesons to two-meson states. They also include the 
exchange or Fock terms for the basic interaction between quarks. 

With a nonlocal interaction between the quarks, there are two-body terms in the symmetry currents ||l8|,|l9|,^ and 
these must be taken into account in the NLO contributions to the couplings of mesons to external currents. We have 
obtained the NLO pieces of the currents by first Fierz transforming our interaction and then applying the technique 
of Refs. ||l^,^. As a consistency check on our treatment we have verified that the GMOR relation holds at NLO. 

In the simple version of the model with only scalar and pseudoscalar interactions we find that the NLO contributions 
to a number of quantities are small for the parameter sets of most interest (the ones with LO quark self-energies of 
about 300-500 MeV). These quantities include the quark condensate and the tadpole piece of the quark self-energy. 
The small size of the NLO pieces is a consequence of cancellations between the effects of virtual mesons and Fock 
terms, as also noted by Ripka p^ . Similar cancellations mean that the NLO contributions to the pion mass and decay 
constant are also small. At NLO these quantities depend not just on the NLO pieces of the quark condensate and 
pion-quark coupling, but also on the tadpole piece of the quark self-energy. The fact that this NLO tadpole energy 
is small is crucial to the smallness of the corrections to pion properties. In other cases, such as the "meson cloud" 
pieces of the quark-self energy, NLO contributions are larger. However they typically contribute at the ^ 25% level 
and so give reason to hope that the 1/Nc expansion is a useful organizing principle for this model. 

In this context, it is worth making a remark about the recent claim of Kleinert and Van den Bosschc ||l^ that 
meson fluctuation effects in the NJL model are large enough to restore manifest chiral symmetry, at least for small 



^This supports a suggestion made in Ref. that the weak coupling of the scalar channel to two pions above the sigma mass 
means that the broad width found for ai an in the model need not be inconsistent with the experimental observation of a 
small width for ai ~* n{nn)a. 
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values of Nc- Although our approach is based on a loop expansion, one might have expected to see large decreases in 
the magnitudes of quantities like the quark condensate, the quark self-energy and the pion decay constant if mesonic 
fluctuations were driving the system towards symmetry restoration. In fact pionic fluctuations do make si gnif icant 
contributions in this direction, but they are largely cancelled by the Fock terms. In the language of Ref. [Q, the 
Fock terms contribute at NLO in 1/Nc to an increased "stiffness" against mesonic fluctuations. This raises a question 
about the estimate of the critical stiffness in that work, which relies on an expression for the stiffness obtained from 
the LO pion propagator. This is in addition to the questions raised by other authors about the universality of the 



critical stiffness in 3 + 1 dimensions |37| and the need to choose an additional cut-off to regulate meson loops in the 
local NJL model |l|,|l§. 

The sigma meson is of particular interest since it can be excited by forces which act to restore manifest chiral 
symmetry, and its mass can be thought of as describing the forces against symmetry restoration. A light sigma seems 
to be a common feature in all NJL- type models at LO. We find that such a state is still present when NLO effects are 
included, although its mass is increased by about 30% (as one might expect from l/Nc arguments). For our preferred 
parameter sets we obtain a sigma mass in range 600-650 MeV. This suggests that the light sigma meson is a real 
feature of this type of model and not just an artefact of LO treatments. 

However, the phenomenological identification of a corresponding low-lying scalar, isoscalar resonance remains con- 
troversial pi|-p3|. As discussed above, the problem is the strong coupling of such a state to the two-pion channel, 
leading to a very broad resonance. In our model we find a width of 100-160 MeV from the LO coupling to two 
pions. While large, this width is much smaller than most phenomenological determinations We have not 

yet calculated NLO contributions to the cttttt coupling, for which terms involving the vector mesons may prove to be 
important. In the extended model at NLO there will be t-channel exchange of a p meson between the pions. The 
attractive force between the pions generated by this mechanism is known to be important for the sigma resonance 
p2| , ^ and could increase its decay width. 

Finally, we have made a preliminary study of NLO effects in the extended model with vector and axial couplings. 
We find that NLO contributions, in particular from the longitudinal vector and axial channels, can improve the high- 
energy behaviour of the quark propagator. In general NLO effects seem to be significantly larger than in the simple 
version of the model, and the convergence of the 1/A^c expansion may be less good as a result. 
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APPENDIX A: PION DECAY CONSTANT AT NLO 

Although the NLO contributions to the pion decay constant arise from a rather large set of diagrams, it is possible 
to deduce some useful cancellations amongst them. A complete description of these cancellations would be somewhat 
lengthy and so we offer only a few brief comments here. Full details can be found in Ref. ||2^ ]. 

Obviously the diagrams that couple the local axial current to the pion all contain an insertion of ^fi'y^- Contracting 
this with g'^ to isolate the longitudinal component one can then substitute for ^75 from the following identity: 

il5 = S-\p+h5 + l5S-\p-) + (m(p+) + m(p_))75. (Al) 

The contributions from the first two terms here can be simplified by cancelling the LO inverse propagators with quark 
propagators that appear in the loop integrals. 

Moreover, for each local diagram there is a similar nonlocal one which contains an additional one-quark loop. (For 
example, compare Fig. ^ with ^ and Fig. Id with |d.) These can be combined with the local contributions, as 
was done at LO ||l^. Each of the NLO nonlocal diagrams can be written as a sum of products of two-quark and 
one-quark loops. In one of these terms, the one-quark loop can be simplified by using the ladder SDE to yield a factor 
proportional to m{0) — rric. There is then a canc ellat ion between this nonlocal term and the part of the corresponding 



in Eq. (Al) 



local di agra m arising from the final term of Eq. ( Al) which leaves only 2mc surviving from the factor m(p+) + m{p-) 



After exploiting Eq. (Al) as outlined above, several other useful cancellations can be identified. As an example of 
the procedure, consider the diagram shown in Fig. When Eq. ( |AlD is substituted for ^75 in the triangular loop 
integral, the first two terms reduce to a sum of two-quark loops. Up to a factor of — GJ the resulting contributions 
are then similar in form to those of the diagram shown in Fig. 
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Hence, in the sum of the contributions from Fig. jija and one can cancel one of the intermediate meson propagators 
to leave a product of a two-quark loop, a three-quark loop and one meson propagator. It turns out that this can then 
be profitably combined with the contribution of Fig. which has the same structure. 

We give here the result of a full analysis of the NLO contributions to f„ in the version of the model with only scalar 
and pseudoscalar interaction s. No te that, rather than defining further loop integrals whose structure is very similar 
to those given above in Sec. [II C , we simply indicate how the products of form factors need to be changed in those 
integrals. The final expression for /tt at NLO is 

fi'Hq')^'^^^4'J{q') I ^Tr[5(fc)]/(fc)(/(fc + ,) + /(fc-,)) 



(0) 

(0) 

(0) 
g-Kqq'fric 

? 



(jW from Eq. (20) with f{v+)f{V-) ^ f{p+)f{P-) 
(jW from Eq. (|2l|) with f{p+)f{p-) ^ f{p+)f{p-) 
(jW from Eq. @ with f{p+)f{p-) ^ f{p+)f{p^) 



in L-^^kn 

^{fip + q) + fip-q) 



§(i-ft,«>(,^))E/|fi|^^«<^-') 



xTr [Sip)T,Sik)T,Sip)] 
>^f{p){f{p + q) + .f{p-q))f{k) 



''-ii fi"\q')- (A2) 

g-nqq 



APPENDIX B: CHIRAL RELATIONS 



A large number of diagrams are required in a NLO treatment of the nonlocal NJL model, as described in Sec. III. It 
is important to have some check that a consistent set has been identified and evaluated. To this end, we demonstrate 
that low-energy chiral constraints are satisfied in the model at NLO. (Again we consider only the simpler version of 
the model, without vector and axial interactions.) 

First consider the pion mass. The chiral expansion of the pion denominator, Eq. (uq), at LO is [18| 



1 - Gi jW(g2) ^ -Gim,M|! ~G,^ + 0{q\m% (Bl) 

where the subscript is used to denote a quantity evaluated in the chiral limit. (We apologize for the fact that the 
notation inevitably becomes rather cumbersome at this point.) If one substitutes the LO chiral expansion into the 

NLO on-shell condition, it is immediately clear that the chiral expansion of jii^(g^) must start at 0{q^, nic) in order 
for the pion to remain a Goldstone boson in the chiral limit. 

A proof of this statement in the local NJL model was presented in Ref. S. Although one has to keep track of 



interaction form factors, the proof for the nonlocal model proceeds along very similar lines |26|. In particular, when 

one takes the chiral limit of the contributions to J^^ from Eqs. ( |2C| ) and (|2^), each of the Dirac traces yields a factor 
that cancels with the denominator of one of the LO quark propagators. The three-quark loop in the contribution of 
the tadpole term to Eq. ( po|) can be reduced to a two-quark integral which can then be used to cancel the denominator 
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of the (t(0) propagator from Eq. (pj|). In the case of the contribution to ji-i"* from Eq. (pT[), the Dirac trace takes an 
even more convenient form in the chiral hmit. For both the sigma and pion exchanges, it factorizes into pieces which 
cancel the denominators from two LO quark propagators. 

An important property noted in the proof of Ref. Q concerns the second diagram of Fig. ^ which involves a virtual 
pion and a virtual sigma meson. In the chiral limit the three-quark loop integrals are proportional to the difference 

between LO polarization loops: Ja^J — . This allows one to replace the product of the virtual meson propagators 

with their difference. It is then possible to combine the contribution from this diagram with other contributions, 

which contain only a single meson propagator. The same property holds in the nonlocal modelR. 

(I) . ... . . 

Many of the simplifications that lead to the vanishing of Ji,/ in the chiral limit can also be exploited in calculating 

its 0{mc) and 0{q'^) terms. One finds that |p6| 



m2(0) ' m2(0) mo(0) 
where the NLO contribution to the condensate is given by: 



(0)3 
9TTqqO 



Substituting the above expansion into the on-shell condition (Eq. (Eq)) yields the pion mass to this order, 



(0) 
qqO 



y-rrqqO 



(mo(0) + Co)2 



(1)1 



0{ml) + 0{N-^). 



(B2) 



(B3) 



(B4) 



We see that the chiral expansion of the pion Bethe-Salpeter amplitude at NLO retains a similar structure to that at 
LO. The changes include the expected ones from the NLO contributions to the quark condensate and the pion-quark 
coupling. However, the form of the dynamical mass scale that appears in Eq. (B4) is less obvious, the shift in this 
being given entirely by the contribution of the tadpole diagram, Eq. ( pj| ) . This is despite the fact that the meson cloud 
of a dressed quark contributes a term of 0{mc) to the chiral expansion of jii^ and also makes a finite contribution to 
the quark self-energy at zero momentum. The fact that the NLO shift in the pion mass is controlled by the coefficient 
of the tadpole diagram has also been noted in the local NJL model P,p^|, in which the overall mass scale can be 
recognized as the stationary point of the one- meson-loop effective action [|6, p!7|. 

From Eq. (84) it can be seen that the Gell-Mann-Oakes-Renner relation p5|| will be satisfied if a version of the 
Goldberger-Treiman relation |38| holds in the chiral limit at NLO, 



f(o) , Ai) 

JnO J-rrO 



(0) 
9TiqqO 

mo(0) 

(0) 
9irqq0 



y-irqqO 
_ gigqO TOo(O) 



_Q 

(0) 
9-nqqO 



(B5) 



18[ to be equal to the LO part of 



The term of order on the right-hand side of this condition was shown in Ref. 

h- ~ 

In order to demonstrate that the NLO part of the decay con stant reduces to the next two terms of Eq. (B5) in the 
chiral limit, first note that all of the contributions in Eq. (A2) are of zeroth order in the chiral expansion. Hence in 
all of the int egra ls the chiral limit may be taken directly, without making an expansion of the integrand. In the first 
term of Eq. ( A2) the ladder SDE can be used to perform the integral, giving 



^It is in fact straightforward to generalize the proof that J^Jq(O) = to the extended version of the model, including other 
mesons. In so doing the products of meson propagators occurring in diagrams of the type shown on the right-hand side of 
Fig. I can be dealt with using the relations jff{q^) - ji°J{q'^) = jj^T (9^) " >^iT(9^) = Jyl^il'^) - >^iT(9^)' where T and 
L denote transverse and longitudinal components respectively. Note, however, that in the extended model, one must take vrai 
mixing into account and so to complete the proof that the pion is a Goldstone boson at NLO, it is also necessary to establish 
that f^l^iO) = 0. 
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■2mo(0)% 



g2 mo(0) g2 ^^(o) ^^(o) 



2Co 



(B6) 



where Eq. (B2) has been used for the chiral expansion of J, 



(1) 



The next three terms in Eq. ( [A2| ) are expHcitly proportional to nic and contain integrals which are very similar to 

those occurring in J^w*, but with different combinations of form factors. Procedures for manipulating such integrals 
in order to simplify their sum were described in Ref. and also above, where they were exploited in order to show 
that jii^ vanishes in the chiral limit. Very similar simplifications can be made here, despite the different form factor 
structures that appear. These give rise to a nonvanishing contribution to /^^ 



(1) 



' 



(0) 
Q-KqqO'^c 

g2mo(0) 



Co /T;\(0) 

(Wo 



mo(0) 



(B7) 



The fifth term in Eq. (A2) is proportional to the tadpole contribution C. Using Eq. (Bl) to expand the factor of 
1 — GiJ^-} {q^), the chiral limit of this term is 



y-n-qqO 



(0) 



mg(0) 



(0)2 
drrqqO , 



(B8) 



In the chiral limit, the integral in the sixth term of Eq. (A2) reduces to the one in the definition of C, Eq. 
giving 



(0) 
d-KqqQ 



Co(l 



r t(o) 



(0) 



(0)2 
5irgg0 / 



(B9) 



The final term of Eq. (A2) is trivially 



fflggO mo(0) 



(BIO) 



Adding together Eqs. (B6) to (BIO) reproduces exactly the 0{N^ ^) terms of Eq. (p35|), demonstrating that the 
Gell-Mann-Oakes-Renner relation holds in our treatment of the nonlocal NJL model at NLO. 
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APPENDIX: TABLES 



mo(0) 




Gi(GoV"2) 


iric 




A 


m(0) 




g-Kqq 


200 




14.3 


4.8 




1459 


245 


246 


2.56 


250 




30.5 


7.8 




1064 


298 


210 


3.13 


300 




53.8 


11.0 




861 


351 


189 


3.70 


350 




85.9 


14.2 




734 


406 


173 


4.28 


400 




128.1 


17.5 




647 


461 


162 


4.87 


450 




181.7 


20.8 




583 


516 


153 


5.46 


500 




248.0 


24.1 




535 


572 


146 


6.04 


TABLE L 


Values of the model parameters, fitted at LO. Also listed are the dynamical quark mass, the quark condensate in 


the chiral limit and the pion-quark couplinj 


g. Apart from Gi 


and gTvqq, which is dimensionless, all quantities are given in MeV. 


mo(0) 






C 




Co 


Co(7r) 




Co (a) 


200 




259 


25.0 




25.3 


90.2 




-64.9 


250 




215 


12.4 




8.0 


44.2 




-36.2 


300 




190 


5.2 




-1.7 


22.4 




-24.1 


350 




174 


0.005 




-8.8 


8.6 




-17.4 


400 




162 


-4.1 




-14.6 


-1.4 




-13.1 


450 




153 


-7.7 




-19.6 


-9.5 




-10.2 


500 




145 


-10.8 




-24.2 


-16.3 




-7.9 


TABLE IL 


Properties of the quark propagator at NLO. All quantities are 


quoted in MeV. The quark condensate is 


evaluated 


in the chiral limit. 


The constant C appears in the part of the self-energy from the NLO tadpole d 


iagram and is 


defined in 


Eq. digl). Its value 


in the chiral limit is Co. 


The contribution to C from the tadpole diagram with a virtual meson m is denoted 


C(m). 


















mo(0) 




(0) 








r(0) (cr tttt) 




Re [mo-] 


200 




385 




1092 




63.5 




455 


250 




423 




1336 




94.4 




552 


300 




454 




1562 




126.3 




624 


350 




477 




1732 




152.0 




650 


400 




492 




1783 




158.7 




>574 


450 




489 




1489 




111.0 




>504 


500 




478 




1515 




116.2 




>450 



TABLE III. Sigma meson mass evaluated at NLO. For comparison the mass and other properties at LO are also given 
(denoted by a superscript (0)). All quantities are quoted in MeV. For the parameter sets with large mo(0), the energy at which 
pseudo-threshold effects are required 13 is rather modest, and hence we quote only a lower bound for the mass at NLO. 
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APPENDIX: FIGURES 




FIG. 1. The LO couplings of a pion to the axial current (represented by the wavy line). 



FIG. 2. Diagrams contributing to the NLO quark self-energy. Double lines are used to denote the LO qq scattering matrix. 



FIG. 3. Meson exchange and two-meson diagrams in the BSE at NLO. There is also a second two-meson diagram (not 
shown) in which the direction of the arrows in one of the quark loops is reversed. 
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(a) 



(b) 




(c) 



(d) 



FIG. 4. NLO couplings of a pion to the axial current involving either meson exchange or two-meson intermediate states. 
There are also two-meson diagrams in which the direction of the arrows in one of the quark loops is reversed. 




FIG. 5. The exchange contribution to the coupling of the pion to the nonlocal axial current. Note that this Fock diagram 
has been distinguished from the similar LO diagram by separating the quark lines to indicate the flow of colour. 





(a) 



(b) 



FIG. 6. NLO contributions to the coupling of a pion to the nonlocal aadal current involving a virtual meson. There is also a 
second diagram similar to (b) in which the virtual meson couples to the antiquark line. 
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FIG. 7. The function a(p) in the quark propagator at NLO, plotted against in GeV^. Also shown are the contributions to 
a(p) obtained by dressing the quark with pion and sigma clouds separately. The parameter set used is the one with mo(0) = 300 
MeV in Table |. 




0.1 - 



-0.25 -0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 

FIG. 8. The function b(p) in the NLO quark propagator, plotted against in GeV^. Also shown are the LO contribution 
and the sum of the LO result and the tadpole contributions. The parameter set is used is the one with mo(0) = 300 MeV in 
Table g 
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FIG. 9. The function a{p) in the quark propagator of the extended model, plotted against in GeV^. The full NLO result 
is shown together with the contributions from the pion and sigma clouds. 
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FIG. 10. The denominator of the scattering matrix in the pion channel at LO and NLO as a function of timelike momentum. 
The parameter set used is the one with mo(0) — 300 MeV in Table ^ 
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FIG. 11. The real part of the denominator of the scattering matrix in the sigma channel at LO and NLO as a function of 
timelike momentum. The parameter set used is the one with mo(0) = 300 MeV in Table |. 
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FIG. 12. The imaginary part (in MeV'^) at NLO of the loop J in the sigma channel as a function of timelike momentum. 
The parameter set used is the one with mo(0) = 300 MeV in Table |. The solid line shows the result of a direct evaluation; the 
dashed line an estimate based on the Cutkosky cutting rules. 
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